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ABSTRACT 






Aims. We present the results obtained with the CoRoT satellite for HD 50870, a S Set star which was observed for 1 14.4 d. The aim 
. ■ of these observations was to evaluate the results obtained for HD 50844, the first 6 Set star monitored with CoRoT, on a longer time 

Ph ' baseline. 

^0 \ Methods. The 307,570 CoRoT datapoints were analysed with different techniques. The photometric observations were complemented 

over 15 nights of high-resolution spectroscopy with HARPS on a baseline of 25 d. These spectra were analysed to study the line profile 
variations and to derive the stellar physical parameters. Some iivhy photometric observations were also obtained to better characterize 
the pulsation modes. 

Q ' Results. HD 50870 proved to be a low-amplitude, long-period spectroscopic binary system seen almost pole-on (/ ^ 21°). The brighter 

5-H , component, which also has the higher rotational velocity (I'sin/ = 37.5 km s"'), is a 5 Set-type variable with a full light amplitude 

0^ ■ variation of about 0.04 mag. There is a dominant axisymmetric mode (17.16 d"'). Moreover, there are two groups of frequencies 

C^ ' (about 19) in the intervals 6-9 and 13-18 d"', with amplitudes ranging from a few mmag to 0.3 mmag. After the detection of about 

250 ternis (corresponding to an amplitude of about 0.045 mmag) a flat plateau appears in the power spectrum in the low-frequency 

region up to about 35 d"' . We were able to detect this plateau only thanks to the short cadence sampling of the CoRoT measurements 

(32 s). The density distribution vs. frequency of the detected frequencies seems rule out the possibility that this plateau is the result 

of a process with a continuum power spectrum. The spacings of the strongest modes suggest a quasi-periodic pattern. We failed to 

find a satisfactory seismic model that simultaneously matches the frequency range, the position in the HR diagram, and the quasi- 

vN , periodic pattern interpreted as a large separation. Nineteen modes were detected spectroscopically from the line profile variations and 

Cn , associated to the photometric ones. Tentative £, m values have been attributed to the modes detected spectroscopically. Prograde as 

>--^ ' well as retrograde modes are present with £ values up to 9. There are no traces of variability induced by solar-like oscillations. 



Key words. Asteroseismology - Stars: variables: 6 Set - Stars: oscillations - Stars: interiors - Stars: binaries: spectroscopic - Stars: 
individual: HD 50870 

1 . Introduction particular, the lowest modes tend to constitute a flat plateau with 

a high-frequency cut-ofF(e.g., around 30 d ' for HD 50844) af- 



,5 Set stars are multimode low-amplitude, short-penod, opacity- ^^^ ^^e simultaneous fit and prewhitening of the strongest ones. 
^ ; driven pulsators in the lower part of the Cepheid instabil- ^j^^ ^ -^^^ interpretation of this fact leads to two main 

- - ity strip. Recent photometric observations performed by the possibilities. This bunch of peaks could be explained by the 

CoRoT (convection. Rotation and planetary Transits) satel- excitation of a large number of modes with a high £ com- 

hte ( |Baghnetal. || 2006|) have shown that about 1000 significant ^-^^^ ^-^^ ^^ insufficient cancellation eff^ect of the flux varia- 

peaks have to be removed to obtain amplitude spectra that show tions when integrated over the whole stellar disk (HD 50844, 



noise only ( |Poretti et al. | 2009t [ Garcia Hernandez et al. || 2009t) . In Ippretti et al.ll2009t iDaszvnska-Daszkiewicz et al.ll2006l) . In this 



respect, iLignieres & GeorgeotI (l2009l) showed how rotation en- 



Send ojfprint requests to: L. Mantegazza ^^^^^^ ^^^ non-cancellation effects, finding that there are three 

e-mail: luciano.mantegazza@brera.inaf.it .. . ., , , . , ^ ^- ^ rv m .u 

. ^ times more visible modes in a star rotating at w = 0.59 w^: than 

The CoRoT space mission was developed and is operated by the . ^ ^. ^ . ^, -i,-!-^ • ^u ^ ^u i ^ 

r^ , ^xTT-o -.u .■ ■ .■ rT-oi' TiooT-f J in a non-rotating star. Another possibility IS that the plateau IS 

French space agency CNES, with participation ot ESAs RSSD and ,, , ,. •■,,, , r 

Science Programmes, Austria, Belgium, Brazil, Germany, and Spain, generated by die granulation noise induced by the surface con- 

This work is based on ground-based observations made with ESO section (|KaUinger & Matthews! |2010D. In connection with the 

telescopes at La Silla Observatory under the ESO Large Programmes possible presence of a convective envelope, |Antoci et al.| (|201l|) 

LP 182.D-0356 and LP 185.D-0056 and on data collected at the reported the detection of solar-l ike oscillations in the 5 Set star 

Observatorio Astronomico Nacional San Pedro Martir (Mexico). HD 187547 observed by JTepfer. iBalona & Dziembowskil(l201 ih 
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Fig. 1. Examples of the hght curve of HD 50870 observed by CoRoT. The lower panel is a zoom-in of a region of the upper panel. The red points 
clustering around zero are the residuals after subtracting the 734 detected frequencies. 



analysed the Kepler data of 6 Set stars and found that in general 
the frequency density is quite moderate. They suggest that the 
high density of the frequencies in the HD 50844 data may be an 
exception. 

Additional detailed analyses of the light variation of S Set 
stars are necessary to clarify these open questions. The aim of 
the present study is add another tile to the assembling mosaic 
by analysing the C oRoT measureme nts of HD 50870 (y=8.88, 
B - y=+0.27, FO; lMcCuskevlll956h. a^ Set discovered in the 
CoRoT preparatory work dPo retti et al.ll20 05\. To achieve the 
most complete picture of HD 50870 we complemented the space 
photometry with ground-based high-resolution spectroscopy and 
with uvby photometry. 



2. The CoRoT data and their frequency content 

CoRoT monitored HD 50870 in the second long run in the di- 
rection of the anticentre (LRa02), started on November 13, 2008 
and finished on March 8, 2009 (AT=1 14.41 d). For the analysis, 
we used the reduced N2 data rebinned at 32 s and only consid- 
ered the 307,570 datapoints for which no problem (i.e., flag=0) 
was reported. The amplitudes of the background orbital varia- 
tions were kept at a very small level by the great effectiveness 
of the baffle, and the subsequent rejection of the uncertain points 
strongly minimized the orbital effects. The light curve was de- 
trended withajow^^degreerolynomial fit to remove the effect of 
ageing ( lAuvergne et al.ll2009h . 

The original data binning of 32 s is higher than necessary for 
the relevant frequency range of 5 Scuti stars and therefore we 
grouped the original data into new bins with a time separation of 
0.005 d to gain CPU time, thus obtaining 20535 datapoints and a 
Nyquist frequency of 100 d'. An example of the resulting light 
curve is shown in Fig.[T] 

The spectral window is shown in Fig. [2] The positions of the 
alias peaks is very similar to those in the window of HD 50844; 
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Fig. 2. Top panel: the spectral window of the CoRoT data of HD 50870. 
Bottom panel: the amplitude spectrum of the original CoRoT data after 
the removal of the slow trend. 



their nature has been carefully discussed bv lPoretti et al.l(l2009l) . 
Note also how the amplitude spectrum reflects the spectral win- 
dow above 40 d"', because of the predominant term at 17. 16 d"' . 
The detection of the frequencies contributing to the light 
variations was pe rformed iterative ly by computing the Fourier 
power spectrum (lDeeming|ll975h . selecting the highest peak. 
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Fig. 3. Distribution of the frequencies vs. amplitudes. Top panel: 
Number of frequencies with an amplitude higher than a given value. 
Bottom panel: histogram of the 734 frequencies (black solid line); dis- 
tribution of the frequencies that were not detected by SigSpec (red 
dashed line). Note the logarithmic scales. 
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Fig. 4. Evolution of the frequency analysis of the CoRoT time se- 
ries. Top panel: local S/N of the n''' highest peak. Middle panel: am- 
plitude of the n* highest peak. Bottom panel: residual rm.s. after 
subtracting all detected frequencies up to n''' term. Note the logarith- 
mic scales on both axes. In the top panel the two points with the largest 
deviations are the Fourier harmonics of the dominant mode. 



adjusting its frequency by a non-linear least-squares fit of all 
previously detected sinusoidal terms plus the new one, comput- 
ing the residuals produced by this fit with respect to the origi- 
nal time series and proceeding to the detection of a new term. 
At each step the amplitude signal-to-noise ratio (S/N) was com- 
puted for the new detected term with respect to the average noise 
in a frequency range of 5 d ' centred around its frequency. The 
iterations were stopped when a local S/N of 3.0 was reached. 
This approach is simil ar to that performed by the Period04 code 
(iLenz & Bregeij|2005l) . but it is made by means of a completely 
automated f77 code, which also allows detecting more frequen- 
cies (Period04 is limited to about 200 components). The fre- 
quency analysis was performed using the relative fluxes, then 
the results were scaled into mmag units. 

We detected 734 frequencies with amplitudes ranging from 
12.961 mmag down to 0.020 mmag. The distribution of the de- 
tected frequencies vs. their amplitudes is shown in Fig. |3] The 
developments of the analysis are plotted in FigH) where we can 
see step by step the amplitudes of the detected terms, their local 
S/N at the detection, and the decrease of the residual rms. These 
plots allow us to see how many components are needed to go 
below a certain level of S/N or rms residual or amplitude. 

Figure |5] shows how the amplitude spectrum changes its 
shape at the various stages of the analysis. At first one strong 
dominant peak is visible at /i = 17.1617 d"\ then two groups 
of peaks show up in the ranges 7-9 d ' and 12-18 d '. 
Successively, bunches of peaks appear between and 35 d"'. 
It is only after the detecting of about 250 terms that the 0-35 d"' 
region becomes a flat plateau, progressively decreasing in height 
after performing new peak detections. The dominant peak at 
about 34 d"' in the third panel from the top is the 2/i harmonic, 
while the isolated peak at about 51 d ' in the fourth and fifth 
panels from the top is the 3/i harmonic. No higher-order har- 



monic of /i was detected. The white noise level is 0.0031 mmag 
as estimated in the region 60-90 d"', which leads to S/N=6.7 for 
the plateau below ~35 d"' visible in the bottom panel of Fig. |5] 
The components with an amplitude higher than 0.3 mmag 
(24 frequencies) are listed in Table [1] Among the detected peaks 
/i5 (13.9725 d"') corresponds to the satellite orbital frequency, 
while /4ig (2.0029 d"', 0.037 mmag) could be produced by the 
gaps introduced in the timeseries by the transit of the satel- 
lite through the South Atlantic Anomaly (SAA), and therefore 
they are probably not attributable to the star. In the solution of 
HD 50870 we find some frequencies with a high rank of detec- 
tion that exactly match the linear combination between a lower- 
rank frequency and the first, predominant term /i (last column 
in Table [1]). They probably arise from nonlinear amplitude dis- 
tortion and are commonly found in multiperiodic radial pulsators 
(Cepheids, RR Lyr, high-amplitude 6 Set stars) as well as in low- 
amplitude 6 Set stars, e.g. FG Vir (Bregeret al. 2005) . The sums 
and, in some cases, the differences are detected for the modes up 
to /lo. The identification of some high-rank frequencies with the 
harmonics of low-amplitude terms (/20, fa, and /24) is reported 
in Table [Has a mere possibility, since they are embedded in the 
flat plateau below 35 d"' and an accidental agreement could oc- 
cur. 

The data were also analysed with the SigSpec code (iReegenI 

l2007h . To go below S/N=4.0 (corresponding to a value of the 
sig parameter of about 5.46), SigSpec requires 790 frequencies. 
The difference in the number of significant peaks with respect to 
the Fourier analysis is probably due to the different approaches 
in computing the white noise level. When comparing the two 
approaches, we found that 673 of the 734 frequencies detected 
with the Fourier analysis are at a distance shorter than 0.0044 d"' 
(half the formal frequency resolution) from the corresponding 
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Table 1. First 24 frequencies identified in the amplitude spec- 
trum (those with an ampHtude >0.3 mmag). The last colum in- 
dicates possible combination terms also present in the data. 



Term 



Frequency 

[d-'] wnz] 



Ampl. Possible combination 

[mmag] terms 



/. 


17.1617 


198.630 


12.961 


.fl6 - 2/l, /365 - 3/l 


/2 


13.0499 


151.040 


2.342 


/87 = /l + fl, /l04 = /l - /2 
/726 = 2/2 


f3 


16.2507 


188.086 


1.736 


M = f\+h 


,/4 


16.0024 


185.212 


1.593 


hi = /l + /4, /l50 = /l - /4 


fi 


8.1366 


94.173 


1.544 


/% = /i + h, fui = 2/5, 

/597 = 3/5 


k 


13.6605 


158.107 


1.313 


/95 = /i + /e, /i90 = fl - /e 


fl 


7.1320 


82.546 


0.974 


/356 =f+ fl, /554 = 2/7, 
/721 = 3/7 


k 


15.1787 


175.678 


0.979 


/l36 = fi+fs 


f9 


17.2805 


200.005 


0.903 


/l67 = /l + /?, /l87 = fg - fl 


fm 


7.6604 


88.661 


0.824 


/lis = 2/10, /l78 = /l + /lO 


fu 


15.0390 


174.061 


0.836 




fn 


16.2473 


188.046 


0.819 




fn 


17.1863 


198.914 


0.692 




/■|4 


15.8344 


183.267 


0.726 




fis 


15.8217 


183.120 


0.707 




fifi 


13.9725 


161.718 


0.792 




fiv 


34.3232 


397.257 


0.501 




flR 


14.0608 


162.740 


0.471 




fi. 


8.7541 


101.320 


0.390 




,/20 


7.8843 


91.253 


0.384 


./ies = 2/20, /284 = 3/20 


fn 


18.4425 


213.453 


0.377 




fii 


8.2945 


96.001 


0.371 




fii 


7.5452 


87.328 


0.324 


fdO = '^fli 


/24 


8.8886 


102.877 


0.308 


fio = 2/24, /505 = 3/24 



frequency detected with SigSpec. In turn, this means that 61 fre- 
quencies detected with Fourier were not detected with SigSpec. 
The first missed detection is the 141-th frequency, whose ampli- 
tude is 0.090 mmag. The histogram with the detections missed 
by SigSpec vs. the Fourier amplitudes is shown in Fig. |3] Since 
the 734 detected frequencies are concentrated in the interval 0- 
35 d', the probability of randomly finding two independent fre- 
quencies in the same resolution bin with the two approaches is 
less than 0.2. 

To check the reliability of the detected frequency terms with 
the highest amplitude, the dataset was subdivided into two equal 
subsets, each with a length of half the original baseline. Each 
subset was analysed independently and the detected frequencies 
were compared with those obtained from the whole dataset. Of 
the first 150 frequencies detected from the whole dataset 146 
were independently detected from the two subsets. The missing 
components are the number 100, 129, 133, and 137. Component 
number 100 has a frequency of 13.0601 d"', and is therefore 
not resolved in the subsets from the second (13.0499 d '). The 
other three missing frequencies were in any case independently 
detected in one of the two su bsets. 

The very recent result bv lAntoci et al.l(l201 ll) concerning the 
detection of solar-like oscillations in the 6 Set star HD 187547 
prompted us to look carefully at the power spectra of HD 50870. 
HD 187547 shows solar-like frequencies in the range of 43- 
73 d"' (i.e., 500-850 //Hz) with amplitudes up to 0.15 mmag, 
while the 6 Set oscillations are confined to frequencies lower 
than about 48 d"' (550 //Hz). The physical parameters of the two 




SO 40 60 80 

Frequency [d"'] 

Fig. 5. Fourier amplitude spectra at different stages of the analysis. The 
flat plateau below about 35 d"' appears after the detection of about 250 
frequencies. The highest peak in the third panel from top is the second 
harmonic of the dominant frequency, while its third harmonic is the 
peak at the highest frequency in 4''' and 5''' panel. 



Stars, as estimated from their spectra and multicolour photome- 
try, are quite similar, which is also supported by the frequency 
ranges covered by the highest-amplitude modes in both stars. We 
should therefore expect to find solar-like frequencies in the same 
region. For HD 50870 all detected frequencies are confined be- 
low 51.5 d"' (~600 //Hz). No significant peak is present above 
this limit. We estimated the white noise level there to be at about 
3 ppm. Only five of the detected frequencies above 35 d ' have 
amplitudes higher than 20 //mag. One of them is the third har- 
monic of the predominant mode and the others are still close to 
the limit of 35 d"' . Moreover, a frequency analysis performed on 
the data with the original binning of 32 s showed that in the fre- 
quency region between 90 and 500 d ' there are no peaks with 
an amplitude higher than 6 //mag. Since there is no regular spac- 
ing, the frequency analysis of the CoRoT data on HD 50870 does 
not supply any hint of the excitation of solar-like oscillations at 
first sight. 

2. 1 . Amplitude stability 

It is rewarding to investigate the stability of the amplitude of 
the predominant mode /i = 17.17 d"'. To do this, we need to 
consider that the modes /y = 17.28 d"' and /i3 = 17.19 d"' are 
close to /i . Figure [T] shows the variation of amplitude recorded 
when the amplitude modulation is calculated with a frequency 
resolution of 0.17 d ', which is insufficient to resolve the three 
components. 

The modulation of the amplitude is clear, but spurious, be- 
cause it results from the beating of the unresolved components. 
We built a synthetic time series using the constant values of am- 
plitudes, frequencies, and phases given in Table [T] for the three 
components and sampled it in the same way as the CoRoT data. 
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Fig. 6. Distribution of the 734 frequencies detected in the CoRoT data 
(top panel) and of the 19 detected in spectroscopy below the pseudo- 
Nyquist frequency of 2 1 .06 d"' (bottom panel). In the top panel the two 
arrows show the second and third harmonics of the predominant mode. 
The insets show a zoom around the predominant frequency. 



Fig. 8. Amplitude behaviour vs. time of a synthetic signal built with the 
frequency terms detected in the region of the dominant mode keeping a 
constant amplitude: there is an excellent agreement with the behaviour 
derived from the observations. 
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Fig. 7. Amplitude behaviour vs. time of the oscillations in the fre- 
quency range [15.5,18.5] d"' (i.e., around the dominant mode). 



The time-frequency analysis of this synthetic signal shows the 
same behaviour of the amplitude (Fig. [8]l, supporting the con- 
clusion that the amplitudes of the /i, /g, and /n modes remain 
constant over the entire time interval covered by the CoRoT ob- 
servations. 



Modes with a lower amplitude were also analysed search- 
ing for the typical random variations of amplitude and phase of 
stochastically excited (solar-like) oscillations, but none showed 
this behaviour. 



3. Ground-based observations 

3. 1 . High-resolution spectroscopy 

The main batch of spectroscopic observations was performed 
during 15 nights between December 15, 2008 and January 8, 
2009 with the HARPS instrument mounted on the 3.6-m ESO 
telescope. La Silla, Chile. The spectrograph was used in the 
high-efficiency EGGS mode, with resolution R ~ 80, 000. We 
obtained 209 spectra in two runs (ten and five nights, respec- 
tively), covering a baseline of 24.3 d, i.e., 4.7 times shorter than 
that of CoRoT. The exposure time was set to 1200 s and the 
S/Ns ranged from 87 to 175 with a median value of 145. The 
spectra were reduced using a semi-autom atic pipeline, written in 
MIDAS and developed bv lRaineil ( l2003h . To study the spectral 
variations induced by pulsation with the highest S/N possible, 
we computed the mean line profiles using the le ast-squares de- 
convolution (LSD) method dPonati et al.lll997l) in the spectral 
region between 4140 and 5615 A, taking care to omit the in- 
tervals containing the //g and Hy lines. The resulting LSD pro- 
files have S/Ns ranging from 1260 to 3900 with a mean value 
of 2406, i.e., the original S/N was increased by a factor of about 
17. In Fig.|9]we show the average LSD profile (top panel, black 
solid line) and the standard deviation of the original LSD profiles 
(bottom panel, black solid line). 
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Fig. 9. Top panel: average LSD mean profiles of the 2008-09 (black Fig. 10. Example of a small region of the observed spectrum of 

solid line) and of the 2010-1 1 spectra (red dash-dotted line). Bottom HD 50870 (top panel) and its separation by the CRES code in the spec- 

panel: standard deviations of the original (solid line) and of the residual tra of the two components (bottom panel), 
(dotted hne) 2008-09 LSD profiles. 



3.2. Evidence for binarity 

The average profile shows the superimposition of two compo- 
nents, one narrow and another broader The close scrutiny of the 
narrow absorption peak did not show any appreciable movement 
with respect to the broader component during the baseline of 
25 d. A non-linear least-squares fit of two Gaussian components 
was performed on the average profile: the centre of the narrow 
component is shifted by about 2.7 km s"' towards the red edge 
of the spectrum with respect to the centre of the broad one. 

To clarify the possible binary nature of this object nine new 
HARPS spectrograms were taken between December 25, 2010 
and January 8, 201 1 and three more between December 18, 201 1 
and January 12, 2012. The LSD profiles were obtained as de- 
scribed above and their average shows that while there is no 
significant difference between the 2010-11 and 2011-12 spec- 
tra, the two components have moved compared to the positions 
in 2008-09 (Fig.|9] top panel). The narrow component is shifted 
by 5.2 km s"' towards the blue edge in the most recent spec- 
tra. There is also a single FEROS spectrogram (/?=48000) taken 
with the FEROS spectrograph attached to the 2.2 m-ESO/MPI 
telescope on January 18, 2003: its LSD profile shows that the 
narrow component is shifted by 0.9 km s"' towards the red edge. 
We can conclude that we observe a spectroscopic bi nary. 

We disentangled the spectra with the CRES code (Illiiicl2004l) 
to obtain a reasonably accurate estimate of the physical parame- 
ters of the two components. To do this we selected some spectral 
regions that contained lines that are sensitive to the temperature 
only or/and to the gravity only in the FEROS spectrogram and in 
the two average HARPS spectrograms. Thes e spectra w e re then 
normalized with the approach described by llUiic et all (l2004t) . 
As an example. Fig. [TOl shows the 2008-09 average spectrum in 
the wavelength range 4487^492 A and the two extracted com- 
ponents. 



Table 2. Physical parameters of the two components 





Start 


Star 2 


Uncertainty 


Teff [K] 


7660 


7077 


±250 


Spectral type 


A8III 


F2II1 




logg 


3.68 


3.74 


±0.25 


[Fe/H] 


0.0 


0.0 


±0.2 


v,„ic [kms"'] 


3.3 


2.8 


±1.0 


V sin i [km s"'] 


37.5 


8.0 


±2.5 



Subsequently, the physica l parameters of the two sta rs were 
estimated with the SME code ( IValenti & Piskunovl 19961) , which 
fits the observed spectrum with a synthetic one derived from 
an interpolation of a grid of models. Estimates of the effec- 
tive temperatures, gravities, metallicities, microturbolent veloc- 
ities, and projected rotational velocities were obtained by av- 
eraging the result of the differ ent spectral region s (Table |2). 
Finally, we used the ATC code ( IStuetz et alj|2006l) to compute 
two synthetic spectrograms with these parameters and combined 
them with the BinMag IDL visualization tool (constructed by 
O. Kochukov) and compared the combined spectrum with the 
observed average spec trograms (Fig. [TTt . The Vienna Atomic 
Line Database VALD dKupka et alfeoOOl : iPiskunov et al.lll995h 
was used throughout. The ratio between the luminosities of the 
two components at 4500 A is 2.0±0.4. These results should 
be considered as preliminary only: additional observations that 
cover more of the orbital period will be able supply firmer re- 
sults. 
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Fig. 11. Example of a comparison between the observed spectrum 
(black dashed line) and the combination of the synthetic spectra of the 
two components obtained with the binmag2 code (red solid line). 



3.3. Line profile variations 

The detection of pulsation modes in the line profile variations 
(LPVs) was performed both with the FAMIAS code (Zima 20()8h 
and w ith the least-squares technique described by iMantegazzal 
(I2OOOI) using the pixel-by-pixel approach. No LPVs were de- 
tected in the component with the narrow lines. On the other hand, 
clear LPVs were detected in the component with the broad lines. 
Figure[T2]shows the least-squares power spectrum obtained from 
the LSD profiles (top panel) and the spectral window of the spec- 
troscopic data. We clearly note that the shorter baseline of the 
spectroscopic observations (24 d) implies a frequency resolution 
worse than that of the CoRoT photometric data. Moreover, there 
are aliases at +1 d"' from the true peak due to the single-site 
observations, and because the observations were performed in 
two separated runs, each spectral peak is flanked by two sym- 
metrical peaks at +0.051 d"' (Fig. [12] inserted panel in the top 
panel). This hampers a comparison of the photometric frequen- 
cies with the spectroscopic ones, since in some cases there are 
several photometric frequencies that might be associated to a sin- 
gle spectroscopic peak. This clearly entails more uncertainty on 
the correct identification of the modes. Moreover, because of the 
spacing between adjacent spectra, the pseudo-Nyquist frequency 
is 21.06 d '. This explains the bunch of peaks around 42 d ' in 
the spectral window (Fig. [12] top panel) and the peaks around 
25 d"' and above 48 d"' in the first power spectrum (middle 
panel). Because of this, attention was paid to avoid misidentifi- 
cations of reflected frequency peaks both above and below this 
value. The comparison of the findings with the two approaches 
allowed us to identify with some confidence 19 oscillation fre- 
quencies (Table [3]l. All frequencies were associated to a photo- 
metric term (first column). There are some ambiguities owing to 
the different frequency resolution of the photometric and spec- 
troscopic data. Where that was the spectroscopic term was asso- 



Fig. 12. Top panel: Spectral window of spectroscopic time-series. The 
high-resolution view of the fine structure of the peaks is shown in the 
insert. Middle panel: Least-squares power spectrum of LPVs of the 
original data. The ordinatae show the percentual data variance reduc- 
tion. Lower panel: Least-squares power spectrum with the 19 detected 
frequencies given as known constituents. The ordinates show the per- 
centual reduction of the residual data variance after the fit with the 
known frequencies. 



ciated to the strongest photometric one that is within the HWHM 
of the centre of the spectroscopic peak. 

The second and third panels of Fig. [12] show the pixel-by- 
pixel least-squares power spectra of LPVs before and after the 
frequency detection of the 19 above quoted terms. Figure [^(bot- 
tom panel, dashed red line) shows the standard deviations across 
the LSD line profile of the residuals obtained by subtracting all 
the detected components. Evidently, many undetected compo- 
nents are still present. 

The amplitude and phase variations of all spectroscopic 
modes run across the whole line profile (from -10 to 90 km s ', 
some examples are shown in Fig.fTTb: this gives us the certainty 
that the pulsational variability has to be ascribed to the fast rota- 
tor. 

The spectroscopic modes are shown in the bottom panel 
of Fig. [6] The pseudo-Nyquist frequency (21.06 d"') must be 
taken into account when comparing them with the photometric 
counterparts (top panel). The amplitudes (the third column of 
Table [3]l are integrated on the whole LSD profile, while their 
S/Ns (fourth column) are computed by FAMIAS. The shape of 
amplitude and phase curves of the dominant mode (17.1617 d"', 
top panel of Fig.fTJb looks quite entangled across the line profile, 
but it seems to be an axisymmetric mode. This probably indi- 
cates the presence of another unresolved mode superimposed on 
the axisymmetric one. Any attempts to fit this peak taking into 
account other nearby peaks (considering adjacent photometric 
detections) did not supply significant improvements because of 
the complexity of the spectral region around /i (see Fig.|6]l. 
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Table 3. Mode detected by the line profile variations listed in 
order of detection. All of them were associated to a photometric 
term. 



Term 



Frequency 



Spectr. 
Ampl. 



S/N 



e,m 



Phot. 
Ampl. 

[mmag] 



/■l 


17.1616 


198.628 


0.245 


14.8 


1,0; 0,0 


12.961 


fu 


15.0390 


174.061 


0.136 


6.5 


4,1 


0.354 


hi 


12.6917 


146.894 


0.136 


9.7 


6,-2 


0.259 


/i% 


16.8485 


195.005 


0.133 


7.3 


6,5 


0.059 


/203 


16.4743 


190.674 


0.129 


6.8 


7,5 


0.064 


,/58 


7.7266 


89.428 


0.113 


8.5 


7,-3 


0.167 


hs. 


17.6884 


204.726 


0.128 


7.9 


6,4 


0.051 


fis 


11.7787 


136.327 


0.117 


8.0 


6,-5 


0.257 


h 


6.5680 


76.018 


0.110 


5.2 


5,-2 


0.125 


k 


15.1787 


175.678 


0.127 


9.0 


4,1 


0.979 


h 


9.4454 


109.321 


0.101 


5.9 


7,-5 


0.125 


f9 


17.2805 


200.005 


0.076 


5.9 


3,1 


0.903 


fu. 


15.3171 


177.280 


0.075 


6.5 


6,0 


0.095 


h9 


8.2442 


95.418 


0.060 


5.2 


9,-5 


0.262 


fM 


8.8886 


102.877 


0.057 


3.0 


8,-1 


0.308 


h 


13.0499 


151.040 


0.073 


5.0 


2,0(1,0) 


2.342 


h. 


10.6011 


122.697 


0.054 


3.4 


7,-3 


0.112 


f2,8 


16.9705 


196.417 


0.058 


4.7 


8,7 


0.046 


/no 


12.5337 


145.065 


0.052 


3.6 


7,1 


0.094 
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Fig. 14. From top to bottom: light, first three moments and A(b-y) 
(comparison minus variable) variations folded in phase according to the 
main pulsation frequency /i = 17. 1617 d"'. 



20 40 60 80 20 40 60 60 



0.005 
0.004 
0.003 
0.002 
0.001 



' — '„ „„„ i^ I I I I I I T i L i I I I I I I I 

ra 0.003 r ' . • 



0.002 

0.001 





I ' '_ .* ' I 



I ' ' ' I ' ' ' I 




20 40 60 80 20 40 60 80 

Doppler Velocity [km s"'] 



Fig. 13. From top to bottom: amplitude (left panel) and phase 
(right panel) variations across the line profile of one axisymmetric 
(17.1617 d"'), two retrograde (7.7262 and 12.6917 d"'), and two pro- 
grade modes (16.848 and 16.874 d"'). 



3.4. Line profile moments 

An independent search for frequencies was performed on the 
timeseries defined by the 0* (equivalent width, EW), P' (ra- 
dial velocity), and 2""* (line width) order moments. For these 
computations the narrow absorption profile (approximated by a 
Gaussian) was subtracted, because the amplitudes of the moment 



variations are dependent on the average profile shape. Because 
the moments are integrated over the whole stellar disc, they are 
more sensitive to low-degree modes. The results are shown in 
Table |4] We see that the first 15 photometric terms were inde- 
pendently detected at least in one of the moment curves. In par- 
ticular, all were detected in the radial velocity variations. 

All modes detected in the zero-order moments were also 
detected in the first-order moment and in general correspond 
to the strongest photometric terms. The only exception is the 
8.8886 d"' mode, which probably suffers from the contamina- 
tion of the 7.8843 d"' mode. That several of the detected modes 
in the first-order moment were not detected in the second-order 
one, and in particular among those with the strongest photomet- 
ric amplitudes, suggests that these could be axisymmetric modes 
(m = 0). The detection of axisymmetric modes is favoured 
for geometrical reasons because the star is seen almost pole- 
on (see paragraph 3.5). We point out /2=13.0499, /3 = 16.2507, 
/4=16.0024, /5=8. 1366, and /6=13.6605d-' as the most promis- 
ing candidates. 

The dominant photometric mode (17.1617 d"') was detected 
in all moment curves. The variation of the first three moments 
phased according to the dominant period are compared with the 
light variations in Fig. [14] However, while /i is the strongest 
term in the zeroth- and first-order moment, it is not the domi- 
nant term in the second-order one, where the dominant mode is 
/8=15.1787 d"'. We cannot exclude that there may be a recip- 
rocal partial contamination caused by aliasing, considering that 
the distance between the two peaks is 1 .983 d"' . Maybe this con- 
tamination could also explain the perturbed shapes of amplitude 
and phase variations across the line profile (Fig. [13] top panels). 
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Table 4. Modes detected by the line moment variations. All were associated to a photometric term. 



Term 


Frequency 


EW 




Rad. vel. 


Line width 


Phot. Amp. 




[d-'] 


U^Hz] 


[10-2 km s-'] 


S/N 


[kms-'] 


S/N 


[kms-']- 


S/N 


[mmag] 


fi. 


6.5680 


76.018 










4.4 


5.8 


0.125 


fi 


7.1320 


82.546 






0.08 


6.3 






0.974 


/26 


7.3304 


84.842 






0.06 


4.4 


3.7 


4.8 


0.280 


/■in 


7.6604 


88.661 






0.04 


3.0 






0.824 


A 


8.1366 


94.173 


0.66 


4.5 


0.12 


9.4 






1.544 


hA 


8.8886 


102.877 


0.47 


3.5 


0.05 


3.4 


3.7 


5.8 


0.308 


h 


13.0499 


151.040 


1.15 


11.4 


0.24 


17.0 






2.342 


ho 


13.3778 


154.835 










3.2 


3.8 


0.260 


h 


13.6605 


158.107 


0.62 


4.7 


0.11 


7.9 






1.313 


fi. 


13.9725 


161.718 






0.05 


3.9 






0.792 


/l7 


14.0608 


162.740 










5.2 


6.2 


0.471 


/■,, 


15.0390 


174.061 






0.04 


3.0 


9.2 


10.8 


0.836 


h 


15.1787 


175.678 


6.46 


6.1 


0.18 


13.7 


12.1 


14.0 


0.979 


fn 


15.8217 


183.120 


0.45 


4.3 


0.07 


5.3 






0.707 


/l3 


15.8344 


183.267 






0.14 


9.7 






0.726 


u 


16.0024 


185.212 


0.55 


5.2 


0.20 


13.5 






1.593 


h 


16.2507 


188.086 


1.04 


9.6 


0.16 


10.9 






1.736 


u 


17.1617 


198.630 


5.01 


43.9 


1.23 


69.1 


5.2 


6.9 


12.961 


fn 


17.1863 


198.914 






0.12 


6.7 






0.692 


f9 


17.2805 


200.005 






0.11 


6.3 


6.0 


8.0 


0.903 



3.5. Mode identification 

It is difficult to accurately fix the stellar physical parameters, 
considering all the uncertainties connected to the disentangle- 
ment of the spectra and the fact that the stellar parallax is un- 
known. However, if we consider the parameters given in Table |2] 
to be sufficiently re l iable, we can derive from the tables by 
IStraizvs & Kurilieng ( Il981h that the primary star (i.e. the 6 Set 
variable) has an A8 III spectral classification and therefore a 
bolometric magnitude of 1.1 ±0.3 mag. Such a star w ould lie 
practi cally in the c entre of the 6 Set instability strip (iBregerl 
120001) . With Eq. 6 in lBregej(l2000l) and using log g and T^ff from 
Table|2]and the above quoted Mboi, we can estimate that the fre- 
quency of the radial fundamental mode (Q=0.033 d) should be 
about 7.4+2.2 d'. 

From the paper by ISchaller et al.l (1 19921) we find, using the 
above derived T^ff and Mboi values, a mass of 2.1+0.3 Mq. 
Combining this with the estimated log^ ~ 3.7, we derive a ra- 
dius of 3.4 ±1.2 Rq, and we adopted this value for the subse- 
quent mode identification. An attempt to derive a more accurate 
value by applying the Baade-Wesselink technique using light, 
colour (b - y variations are shown at the bottom of Fig. [141 see 
Sect. 3.6 for details) and radial velocity variations of the domi- 
nant mode (17.17 d"') did not lead to a significant result. In turn, 
this supports the ^ = 1 identification for the /i mode. Only for 
this £ value the projected area va riation is null and the Baad e- 
Wesselink technique breaks down (IStamford & Watsonlll98lh . 

The mode i dentification was performed with FAMIAS in the 
AP mode (see IZimal [20081) . Since the profiles shapes are com- 
plex because of binarity, their average value was subtracted from 
the individual LSD profiles and the identification was performed 
on these residuals. A first identification was performed individ- 
ually on each detected mode by fitting its amplitude and phase 
variations across the line profile. This supplied preliminary {,m,i 
values. After that we performed a simultaneous fit of all modes 
with the previously detected £,m values kept fixed, but leaving 
the inclination as free parameter. From the values of the discrim- 



inant we were able to estimate / = 21° ± 6° as the most probable 
inclination of the rotational axis (Fig. [TsT l. Finally, a separate 
fit was performed on each detected mode assuming this inclina- 
tion value to improve the (!,m values. The derived (!,m values are 
listed in Table[3] where the negative m values indicate retrograde 
modes. Uncertainties are estimated to be ±1 for the { degree, 
while they are larger for the m values. The main cause for un- 
certainties are the contaminations by adjacent unresolved modes 
and by aliases. But the distinction between prograde and retro- 
grade modes through the shape of the phase curves is quite un- 
ambiguous (see Fig. [T3T l and the total amplitude of these curves 
across the line profile supp lies a good constraint to the £ value 
(iTelting & Schriiverslll997l) . 

Knowing the inclination, we were able to estimate Veg = 
116 ± 38 km s"' from the vsin/ value and consequently the 
rotational frequency, Vrc = 0.67 ± 0.32 d'. We applied the 
same approach to the secondary component and estimated that it 
could be an F2 III star, with a bolometric absolute magnitude of 
1.8 ±0.3 mag, located just outside the red border of the instability 
strip. 



3.6. Multicolour piiotometry 

We also observed HD 50870 in Stromgren uvby photometry to 
have an independent tool to propose at least an identification of 
the modes with the highest amplitudes, uvby photometry was 
performed at San Pedro Martir Observatory from 2007 to 2010 
and the final dataset consists of 593 datapoints collected on 28 
nights for a total survey of 156.5 h. The error on a single mea- 
surement is 0.005 mag in the vby filters and around 0.015 mag 
in u filter. The frequency analysis of this dataset allowed us to 
detect the terms from /i to /s in an unambiguous way. The 
aliasing effect is noticeable after the sixth term, since instead 
of /6=13.66 d"' the highest peak is the alias at /e+l d"'. The 
subsequent residual spectrum became very noisy and the highest 
peaks cannot be identified in a reliable way with the frequencies 
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Fig. 15. Behaviour of the discriminant obtained by fitting all 19 spec- 
troscopic modes simultaneously to the LSD profiles. 
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Fig. 16. Histogram of the differences between couples of frequencies 
with amplitudes higher than 0.1 mmag in the 0-30 d"' interval. The 
pattern at 7.5, 15, and 22 d"' is indicated. 



detected in the CoRoT timeseries. Looking at the photometric 
amplitudes listed in Table 3, we infer that we were able to de- 
tect from ground all terms with a CoRoT amplitude higher than 
1.0 mmag. 

We calculated the amplitude ratios and phase shifts by us- 
ing th e y colour as reference system (Ganido 2000; Pore tti et al.l 
l2009l) . The predominant term f\ is characterized by negative 
phase shifts (e.g., 0,, - (f>^ - -0.05 + 0.03 rad, (^^ - (f)^, - 
-0.04 + 0.03 rad). The amplitude ratios are A^jAy = 1 .39 + 0.04 
and AblAy = 1.19 + 0.05. Therefore, the case of HD 50870 is 
different from that of HD 50844, where the positive phase shifts 
su pported the identifi cation of f\ as radial mode (see Fig. 10 
in iPoretti et al.ll2009l) . Here, the negative phase shifts suggest a 
nonradial mode, but this identification is ambiguous since the f\ 
mode might belong to a high radial order (we obtain Q=0.014 d 
from the procedure described in Sect. 3.5) and the th eoretical 
phase shifts are negative in this case, too (see Fig. 13 in iGarridol 
I2OOOI) . The phase shifts of the terms from fi to /& also have neg- 
ative values, but the error bars are at least three times larger than 
that the f\ ones, which precludes a reliable mode identification. 

4. Frequency spacings 

The asymptotic regime of the solar- like oscillations describes the 
well-known periodic pattern characterized by the large 



and the small 



5v ■ 



Vn,e - Vn-l/+2- 



(1) 



(2) 



separations. This regime is not valid for 6 Set stars and we 
were unable to detect the two separations as sharp structures 
in the frequency distributions; instead they ap peared as quasi- 
regular spacings (iGarcia Hernandez et al.ll2()()9l) . These spacings 
become more evident when the search is restricted to the modes 



covering a small ( range. Because the modes with C < 5 should 
show the highest amplitudes, we performed our analysis by con- 
sidering the frequencies with amplitudes higher than 0. 1 mmag. 
Searching for the large spacing (Fig. [16) we found a promi- 
nent peak around 7.5 d"' and other peaks at 15 d"' and 22 d"', 
i.e., about double and triple frequency of the former value. This 
suggests a regular pattern with a frequency separation of 7.5 d"' . 
We also computed the Fourier power spectrum of a series of 
Dirac's 6 functions with the same amplitude and centred on the 
frequencies detected in the CoRoT data. This power spectrum is 
characterized by a broad peak centred at about 7.8 d"' (Fig.fTTb. 

When searching for a smaller regular spacing in the 0-5 d ' 
interval, we were able to identify a prominent peak at about 
1.2 d"' (Fig. [T8l l. The 1.2 d"' value is approximately twice 
the rotational frequency (0.67 d ') and corresponds to the peak 
shown in Fig. [17] The recurrency of the 0.6 d ' separation is con- 
firmed by some peaks reproducing a multiplet pattern, as those 
observed at 6.70, 7.13, 7.66, 8.14, and 8.76 d"'. 

The study of the frequency spacings provided two main 
results. First, the Fourier transform and the histogram in the 
0-5 d"' range (Fig. [TTj and Fig. [TST l show the double of the 
rotational splitting, which is predicted for a frequency range 
hig h enough to render th e eff'ect of the Coriolis force negligi- 
ble (iLignieres et al.ll2010l) . Second, the frequency differences in 
the 5-30 d~' range show a s pacing periodicity (as defined by 
iGarcia Hernandez et al]l2009h of about 7.5 d'. 



5. Asteroseismic modelling 

The availability of the physical parameters and of the full 
description of the pulsational spectrum allowed us to per- 
form the computation of asteroseismic models (i.e., equilib- 
rium plus pulsation models), the analysis of instability, and a 
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Fig. 17. Power spectrum derived from the distribution of the frequen- 
cies with amplitudes higher than 0. 1 mmag. The strongest peak is close 
to the first peak of the pattern shown in Fig. [16] The peak at 1.15 d"', 
close to twice the rotational frequency, is also indicated. 
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Fig. 18. Histogram of the differences between couples of frequencies 
with amplitudes higher than 0.1 mmag in the 0-5 d"' interval. The peak 
at 1.2 d"' and its double are visible. 



study of quasi-periodicities in the observed oscillation spectrum 
dGarcia Hernandez et"aDl2009l) . 



5. 1 . The asteroseismic models 

Equili brium models were computed with the evolutionary code 
CESAM (lMorellll997tlMorel & Lebretonll2008l) . Physics included 
in the models was appropriate to the description of 6 Set stars 
(ICasas et al.l l2006|) . We considered a variation in the convec- 
tive efficiency parameter umlt up to 1.5, of the overshoot up 
to 0.3, and a solar metallicity. Diffiision and ra diative forces 
were assumed to be negligible for this type of stars (iGoupil et al.l 
|2005|) and were not included in the modelling. We consid- 
ered about 2000 points for the computation mesh, following 
the studies of the Evolution and Asteroseismic Tools Ac tivities 
(ESTA/CoRoTQ; iMova et al.ll2008l: iLebreton et al.ll2008h . 

We considered / = 21°+6° and then Veq = 116+38 kms"', as 
resulting from the spectroscopic mode identification. For these 
rotational velocities it is expected to find a non-rigid rotation in 
the stellar interior, in particular at the interface between the con- 
vective core and the radiative envelope. In that region the mixed 
modes (low-order p and g modes) have a si gnificant sensitivity 
to variations of the internal rotation profile (ISuarez et al .1120091) . 
Without a detailed description of the shape of the internal rota- 
tion profile, we can reasonably approximate it with the physical 
assumption of a loca l conservation of t he angular momentum 
during the evolution (ISuarez et al.ll2006h in the radiative zones, 
and a quasi-instantaneous transport of the angular momentum in 
the convective zones. The total angular momentum is conserved 
(no mass loss is expected during the main-sequence stage). 

Pulsational models were obtained by a perturbation of the 
equilibrium m odels using the adiabatic oscillation code fslov 
(ISuarezll2002l and references therein). Models were computed 
using the physical parameters obtained from HARPS spec- 
troscopy (see Table|2]i. The uncertainties in those parameters de- 
fined our space of validity for the representative models. Models 
were constrained to predict the highest amplitude modes having 
the expected spherical degrees (Table|3]l. This procedure defined 
two sets of models with masses of M=2.10 and 2.I8M0. 



5.2. The instability analysis 

A non-adiabatic analysis of the mode energy bal ance of the 
sets of models was performed with the GraCo code (iMova et al.l 
120041: iMova & Gairidol |2008|) . Theoretical results were com- 
pared with the range covered by the spectroscopic modes and 
by most of the highest amplitude photometric modes, i.e., 
[6.5,18.4] d ' (Tables [1] and [3]l. Each theoretical model is rep- 
resented by a point in Fig. [19] We found a few models that 
predicted the over-stable region [6.5,18.4] d"'. The analysis of 
instability thus constrains the representative models to a first 
set of models with M=2.18Mq, 3.8692< log T^ff <3.8811, 
3.7003< logg <3.7746, 1.5009< L/Lq <1.5056 and a sec- 
ond set with M=2.10Mo, 3.8673< log T^b <3.8803, 3.7935< 
log g <3.8529, 1.3812< L/Lq <1.3886. 



5.3. Analysis of quasi-periodic patterns 

Using the aforementioned sets of theoretical models, we looked 
for quasi-periodic patterns. There is no model with a mean Av 
around 7.8 d"' (Fig.fTTI). The pulsational models taking into ac- 
count the rotational effects up to the second order of the pertur- 
bative approach did not point out structures similar to the multi- 
plet centred at 7.66 d ' (Sect. |4]i. Consequently, we considered 
the half value of 3.9 d '. We estimated the uncertainty of this 
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Fig. 19. Bottom versus upper limit of the instabihty range. Models 
computed using the physical parameters from spectroscopy are plotted 
as red and green crosses. The limits of the frequencies with the high- 
est amplitudes observed in HD 50870 are depicted with a blue star. The 
shaded yellow area shows where a theoretical model must be to predict 
an over-stable frequency range containing the observed one. The mod- 
els with a Av in the observed range of [3.5,4.5] d"' are marked as light 
blue squares. 



value using the FWHM of the highest peak of the Fourier trans- 
form (Fig. [TTI i. Thus, we obtained a [3.5,4.5] d"' range for Av. 
IVlodels with 2.10 Mq show these values, but they also predict a 
lower limit for the excited modes that is much higher than the 
observed one (Fig.[T9]l. 

Therefore, there is no model that at the same time pre- 
dicts the [6.5,18.4] d"' range of observed frequencies and the 
[3.5,4.5] d"' range of the quasi-periodicity intepreted as a large 
separation. This seems a general problem in the era of rich pul- 
sational spectra detected b y space observations, as stressed by 
lUvtterhoeven et all ( 1201 ll) when analysing the 6 Set stars ob- 
served with Kepler. 

6. Discussion and conclusions 

The space-based photometric and the ground-based spectro- 
scopic monitoring s of another 6 Set star after HD 50844 
dPoretti et al.ll2009h confirm that these variables show a compli- 
cated amplitude spectrum. HD 50870 has one predominant pul- 
sation mode (17.1617 d"', 12.96 mmag), and about 20 modes 
with amplitudes beween 2.3 and 0.3 mmag and frequencies be- 
tween 7 and 19 d'. Well-isolated peaks were detected down to 
0.05 mmag. The frequencies with amplitudes lower than about 
50 /umag (S/N~ 12) form a flat low-frequency plateau with a cut- 
off at about 35d"' (~400 /iHz). After the detection of several 
hundreds of frequencies (down to amplitudes of about 12 ppm), 
there is still an excess of power at low frequencies (see last panel 
of Fig. 4). About 1800 frequency terms are necessary to obtain a 
flat power spectrum by means of a classical Fourier solution of 
the light curve. 
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Fig. 20. Frequency analysis of the HD 50870 light curve but sampled 
with the data binning of the Kepler satellite. Only 145 peaks are de- 
tectable above the S/N level of 4.0. The spectral window shifted at the 
frequency of the predominant mode (17.16 d ') is shown in the top 
panel. 



The results obtained for HD 50870 confirm that we need 
hundreds of frequencies to reproduce the light curves of 6 Set 
stars, i.e., that HD 50844 is no t an exception or a unique case 
(iBalona & Dziembowskill201 lb . Because the light curves of 6 
Set stars show rapid variations, the detailed investigation of the 
signal content must be performed on the basis of time series with 
a short cadence sampling. To illustrate the point, we grouped the 
CoRoT measurements of HD 50870 in bins of 30 min, i.e., the 
typical long cadence of the Kepler data. The resulting ampli- 
tude spectrum obtained from the frequency analysis of the 5135 
datapoints is shown in Fig. |20] The spectrum is generally flat 
and a limited set of frequencies (145) is above the significance 
limit of 4.0. This can be explained by the fact that the pseudo- 
Nyquist frequency of the almost equally-spaced Kepler-like 
data is about 23 d"'. Consequently, a consistent fraction of the 
spectral power below this limit (where most of the stellar signal 
is concentrated) is reflected above it, thus simulating a higher 
noise level. Indeed, we calculated a white noise level of 15 yumag 
(to be compared to the 3 yumag level obtained from the analysis 
with the CoRoT data sampling). We can conclude that the ca- 
dence of 30 min is unable to describe the light variation of 6 Set 
stars in a complete way. 

We remain with the still unsolved problem whether the ex- 
cess of power detected in the CoRoT timeseries below 50 d ' is 
due to the excitation of a huge number of modes or is the result of 
a broad-band low-frequency noise process, i.e., we are observ- 
ing a discrete or conti nuum power s pectru m. A similar plateau 
has been observed by iPoretti et aP (l20()9h in HD 50844 and it 
seems to be present in HD 174936, too (iGarcia Hernandez et al.l 
l2009h . The analysis of the light variations of other 6 Set stars 
observed by CoRoT seems to support this result: HD 170699, 
for instance, which is a very fast rotator (v sin i-27Q km s"') has 
a plateau with the cut-off at about 25 d"' (Mantegazza et al.. 
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Fig. 21. Density distribution of the 1800 frequencies necessary to 
whiten the power spectrum (number of frequencies per d"') as a func- 
tion of frequency. 



in preparation). A tentative preliminary answer can be given for 
HD 50870. The histogram that shows the distribution with a bin- 
ning of 1 d ' of the frequency terms needed to whiten the power 
spectrum (i.e., to lower its mean value at the low frequencies 
at about the same level as above 60 d"') is shown in Fig. |2T| 
Excluding the peak close to zero frequency, which might be af- 
fected by the residual satellite low-frequency noise, the maxi- 
mum density is 63 frequencies per d"', with an average value 
below ~ 35 d ' of about 50 frequencies per d" ' . It is well known 
that if we have a time series generated by a process with a con- 
tinuum power spectrum, and we estimate its power spectrum by 
means of the discrete Fourier transform, we approximate it with 
a discrete spectrum consisting of a bunch of peaks with a separa- 
tion between t hem of the orde r of the reciprocal of the time series 
baseline (e.g.. lChatfieldll984l) . The baseline of the CoRoT obser- 
vations is about 1 14 d, then we should expect a frequency density 
of the order of 1 14 frequencies per d ' to be detected if there is 
a signal with a continuum power spectrum. Since the observed 
density is decidedly lower, this analysis rules out the possibil- 
ity that the plateau is due to a process with a continuum power 
spectrum. However, it does not rule out the action of the gran- 
ulation effect, since this effect could be described as a random 
process with characteristic time scales that could generate dif- 
ferent densities. The observed frequency ranges, amplitudes, and 
densities pose observational constraints to the models of the con- 
vective layers in S Set stars, at the moment limite d to tempera- 
tures cooler than that of HD 50870 (T^g <6850K, ISamadietal] 
I2OO2I) . A more detailed theoretical investigation is necessary to 
support the effectiveness of convection in 6 Set stars. 

As for HD 50844, high-resolution spectroscopy was of 
paramount importance to improve the physical scenario. The 
analysis of the line profile variations allowed us to identify 
nineteen progrades, axisymmetric and retrogrades modes hav- 
ing £ <9 as well as an inclination angle of about 21°. The corre- 



sponding rotational velocity is about 116 km s"', and therefore 
HD 50870 is rotating at about 34% of the break-up velocity. 

HARPS spectroscopy suggests that the predominant mode 
is axisymmetric (m-O), but it does not allow us to distinguish 
between a radial or a nonradial mode. The greater difference be- 
tween a radial mode {{ - 0) and an axisymmetric nonradial one 
with £ = 1 is in the equatorial region, but in the current case 
of a star seen almost pole-on this region supplies a very limited 
contribution through the Doppler effect to the overall line profile 
variations. On the other hand, the Baade-Wesselink technique 
(Sect. 3.5) and the amplitude ratios and phase shifts in differ- 
ent photometric colours (Sect. 3.6) seem to indicate a nonradial 
mode. More in general, despite the huge observational effort, 
there is still a mismatch between the observational contraints 
and the theoretical modelling. The mechanism of selecting the 
modes with the highest amplitudes is not yet fully understood, 
als o when considering t he effect of fast or moderate rotation. 

lAntoci et al.l (l201lh observed a 0.1-0.2 mmag level of the 
solar-like oscillations for HD 187547. We did not find any ev- 
idence of the presence of oscillations like this in HD 50870, at 
least above the level of 12 ppm. The non-detection may be the 
result of the dilution of the structures broadened by the mode 
lifetime into the noise and an apparent reduced amplitude due 
to the light from the companion. In both cases, the damping 
should be very effective to produce a null detection. HD 187547 
has physical parameters very similar to those of HD 50870, but 
also a slight overabundance of some chemical elements. Since 
the amplitude of solar-like oscillations depends on the surface 
metal abundance dSamadi et al.ll2010al lbh. this particularity can 
explain why the visibility of these oscillations is enhanced in 
HD 187547. A well-hidden cooler companion could also be an 
explanation for the case of HD 187547 because we were able to 
detect the binarity of HD 50870 only thanks to the continuous 
monitoring with high-S/N HARPS spectra. 
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